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First-principles study on physical properties of a single ZnO monolayer with 

graphene-like structure 
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The elastic, piezoelectric, electronic, and optical properties of a single ZnO monolayer (SZOML) 
with graphene-like structure are investigated from the first-principles calculations. The phonon 
dispersion curves contain three acoustic and three optical branches. At F point, the out-of-plane 
acoustic mode has an asymptotic behavior Lo{q) — B(f with B — 1.385 x lO"'^ m^/s, while two 
in-plane acoustic modes have sound velocities 2.801 km/s and 8.095 km/s; the other three optical 
modes have frequencies 250 cm~^, 566 cm~^, and 631 cm~^. The elastic and piezoelectric constants 
are obtained from the relaxed ion model. It is found that the SZOML is much softer than graphene, 
while it is a piezoelectric material. The electronic band gap is 3.576 eV, which implies that the 
SZOML is a wide band gap semiconductor. Many peaks exist in the linear optical spectra, where 
the first peak at 3.58 eV corresponds to the band gap of SZOML. 



I. INTRODUCTION 

Many ZnO nanostructures such as nanowires, 
nanobelts, nanorings and nanohelixesi"— have been 
synthesized in recent years, which are expected to 
be used to make novel electronic devices such as 
field-effect transistors, nanogenerators and so on4i^ 
Vievifcd from the small scales, these nanostructures are 
locally wurtzite type. In 2005, Claeyssens et al^ found 
that the graphite-like structure was more energetically 
favorable than the wurtzite structure for a very thin 
ZnO film through the density functional theory (DFT) 
calculations. They obtained the critical number of ZnO 
monolayers in the thin film to be 18, beyond which the 
graphite-like structure would transit to the bulk wurtzite 
structure^. Subsequently, the present author and Hu^ 
pointed out that DFT calculations might overestimate 
the interlayer interactions because the widely used DFT 
packages cannot accurately describe the non-covalent 
bonds^. The critical number was estimated to be 3 
based on some arguments for the interlayer interactions 
between two ZnO monolayers and the DFT calculations 
for a single ZnO monolayer (SZOML) ^ In 2007, Tusche 
et alM- observed graphite-like ZnO monolayers in the 
experiment and found that critical number to be 3 or 4. 

We expect that this new material might have good 
physical properties. However, there is very few investi- 
gation on this topic except our previous wor]^ on elas- 
tic constants and the workiSiii on metallic edges and 
magnetic behavior in ZnO nanoribbions. Here, we will 
address the elastic, piezoelectric, electronic, and opti- 
cal properties of an SZOML with graphene-like structure 
by first-principles calculations (the ABINIT package^^). 
The rest of this paper is organized as follows: In section 
im we show an optimized structure of the SZOML. In 
section IIII[ we calculate the phonon dispersion relation. 
In section lTVl we obtain the elastic and piezoelectric con- 
stants based on the relaxed ion model. In section |Vl the 
electronic band structure and the band gap are obtained 
through the DFT calculation and GW approximation. 
In section I VII we show the results of the optical dielec- 



tric functions. We summarize the main results and infer 
the possible physical properties of single- walled ZnO nan- 
otubes based on these main results in section IVlIl 



II. STRUCTURE 

The DFT calculations^iii^iii and the experiment^ re- 
veal that an SZOML with graphene-like structure is 
chemically stable. Its structure is schematically depicted 
in Fig. [1] where ai and a2 are the two lattice vectors. 
Each unit cell contains one Zn atom and one O atom. 
We optimize the structure by taking Troullier-Martins 
pseudopotentials,^^ plane- wave energy cutoff 50 Hartree, 
and 8x8x1 Monkhorst-Pack k-points^^ in Brillouin- 
zone. The exchange-correlation energy are treated within 
the local-density approximation in the Ceperley- Alder 
form^^ with the Perdew-Wang parametrization.-'^^ We 
keep the lattice constant in z-direction 40 bohr such that 
our result is available for the SZOML. We obtain the two 
lattice constants in xy-plane |ai| = |a2| = 6.062 bohr 
(=3.208 A), and correspondingly, the bond length 1.853 
A, which is close to the previous value^ obtained by the 
present author and Hu with different ecut and k-points. 



III. PHONON DISPERSION RELATION 

The phonon dispersion relation uj = w(g) entirely re- 
fiects the lattice dynamic behavior of SZOML, where lu 
and q are the phonon frequency and wave number, re- 
spectively. Here, we obtain the phonon dispersion re- 
lation through the first-principles calculations following 
the worki2i2£ by Gonze and others and show the results 
in Fig. H 

Three phonon dispersion branches originating from F 
point of the Brillouin zone correspond to the acoustic 
modes: an out-of-plane mode (the lower branch), an in- 
plane tangential mode (the middle branch) and an in- 
plane radial mode (the upper branch) . Among them, the 
phonon dispersion curve of the out-of-plane mode near 
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them is 

kc = (tB'', (1) 

where a — 1.516 x 10^^ kg/m^ is the density (i.e., mass 
per area) of SZOML. Then we obtain = 0.182 eV, 
which is much smaher than the bending rigidity (1.62 
eV}^ of graphene. 

The remaining three branches are optical modes: one 
out-of-plane mode and two in-plane modes. At F point, 
their frequencies are non- vanishing, and the correspond- 
ing values are 250 cm^^, 566 cm^^, 631 cm~^, respec- 
tively. 



FIG. 1: (Color online) Graphene-like ZnO monolayer, ai and 
32 are the two lattice vectors. 
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FIG. 2: (Color online) Phonon dispersion curves plotted along 
high symmetry directions in the Brillouin zone. 



r point has an asymptotic behavior ui{q) = Bq^ with 
B = 1.385 X 10"'' m^/s. The in-plane tangential and ra- 
dial modes have linear behaviors near F point. The sound 
velocities are — 2.801 km/s and Vr — 8.095 km/s, re- 
spectively, which are much smaller than the correspond- 
ing velocities 15 km/s and 24 km/s for graphenci^ Here 
we should emphasize that one cannot directly derive the 
in-plane rigidity and shear modulus from the sound ve- 
locities Vt and Vr because the elasticity and piezoelec- 
tricity of SZOML couples with each other for in-plane 
deformations.— But the piezoelectricity and the out-of- 
plane mode are decoupled. Thus we can derive the bend- 
ing rigidity kc of SZOML from B. The relation between 



IV. ELASTICITY AND PIEZOELECTRICITY 

We have mentioned that the elastic constants of 
SZOML cannot be directly derived from the phonon dis- 
persion relation because the coupling between elasticity 
and piezoelectricity. However, we can directly calcu- 
late the elastic and piezoelectric constants by using the 
ABINIT package. Because the SZOML has threefold ro- 
tation symmetry around 2;-axis and a reflection symmetry 
with respect to the y-axis as shown in Fig. [U the elastic- 
ity and piezoelectricity for zero external electric field can 
be expressed in the matrix form— as 
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respectively. In the above expressions, ui and (T2 are the 
normal stresses along x and y directions, respectively, si 
and S2 are the normal strains along x and y directions, 
respectively, gq and se are the in-plane shear stress and 
shear strain, respectively. is the spontaneous polar- 
ization in the y-direction. Pi and P2 are polarizations 
along X and y directions, respectively, cn and C12 are 
two independent elastic constants, while di is the inde- 
pendent piezoelectric constant. Because the thickness of 
SZOML is unknown, we let it be implicit in cn, C12, 
and di. Based on the relaxed ion modelj^S we obtain 
Cll = 56.8 eV/ZnO, C12 = 40.3 eV/ZnO, and di = 3.0 
Pe/ZnO from DFT calculations, where pe — 8.5 x 10^"^° 
C m is the atomic unit of electric dipole moment. 

Using the elastic constants, we can derive the in-plane 
Young's modulus (Y), shear modulus (G), and Poisson 
ratio (ly) of SZOML as follows: 

Y = cii(l - cl^/cfi) = 28.1 eV/ZnO, (4) 
G = (cii-ci2)/2 = 8.2eV/ZnO, (5) 
= C12/C11 = 0.71. (6) 
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The Young's modulus and shear modulus of SZOML are 
much smaller than those of graphene, 115.4 and 49.6 
eV/CC, respectively, while the Poisson ratio of SZOML 
is much larger than that of graphene, O.lGi^^^ 



V. ELECTRONIC BAND STRUCTURE 

As a new material, we expect that SZOML has a 
good electronic property. Its band structure is calculated 
within the DFT framework. In Fig. [31 we plot energy dis- 
persion curves for 9 valence and 5 conduction bands along 
high symmetry directions in the Brillouin zone. There 
are 1 low energy valence band, 5 middle energy valence 
bands, and 3 high energy bands among the 9 valence 
bands. The direct gap between the lowest conduction 
band and the highest valence band is Egap — 1.762 eV. 
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shift" 1.813 eV (i.e., the correction of band gap obtained 
from GW approximation). 
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FIG. 4: Dielectric functions. 'Im' and 'Re' represent the imag- 
inary part and real part, respectively, h is the unknown thick- 
ness of SZOML. 

In Fig. |4l we show the frequency dependent real parts 
and imaginary parts of dielectric functions eyy{u!) and 
ezz(w), where h is the unknown thickness of SZOML. exx 
and Cxy are not shown because their shapes are similar to 
eyy. The peaks in the linear optical spectra lm{heyy} — u 
exist at 3.58 eV, 3.80 eV, 3.93 eV, 4.20 eV, 4.38 eV, 
4.52 eV, 4.87 eV, and so on, which can be identified 
from the band structure. In particular, the first peak 
at 3.58 eV marked with an arrow in Im{/iej,j,} — w curve 
corresponds to the band gap, 3.576 eV, obtained from 
GW approximation. 
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FIG. 3: (Color online) Energy dispersion curves for 9 valence 
and 5 conduction bands plotted along high symmetry direc- 
tions in the Brillouin zone. 

The large underestimation of the band gap is a well- 
known problem in DFT calculations. We use the 
GW approximation2Si21 to obtain the correct band gap 
^gap = 3.576 eV which is larger than the calculated band 
gap (2.6 eV)2^ of the bulk ZnO. Thus the SZOML can 
be regarded as a wide band gap semiconductor. 



VI. OPTICAL DIELECTRIC FUNCTIONS 

The frequency dependent optical dielectric functions 
of semiconductor can be computed from the first 
principles.—— In terms of the symmetry of SZOML, 
the non- vanishing dielectric functions are txx^ ^yy^ ^zz, 
^xy = ^yx, which can be calculated by taking 81 bands, 
19 X 19 X 1 Monkhorst-Pack k-points, and an "scissor 



VII. CONCLUSION AND DISCUSSION 

We have obtained the phonon dispersion relation, the 
elastic and piezoelectric constants, the electronic band 
structure, and the optical response functions from the 
first-principles calculations. The main results are listed 
as follows. 

(i) The phonon dispersion curves contain three acous- 
tic and three optical branches as in Fig [2l At F point, 
the out-of-plane acoustic mode has an asymptotic be- 
havior uj{q) = Bq^ with B = 1.385 x 10"'' m^/s, while 
two in-plane acoustic modes have sound velocities 2.801 
km/s and 8.095 km/s; the other three optical modes have 
frequencies 250 cm^^, 566 cm~^, and 631 cm~^. 

(ii) The elastic and piezoelectric constants are cu — 
56.8 eV/ZnO, Ci2 = 40.3 eV/ZnO, and di = 3.0 Pe/ZnO. 
The SZOML is a piezoelectric material. The in-plane 
Young's modulus, shear modulus, and Poisson ratio are 
calculated as 28.1 eV/ZnO, 8.2 eV/ZnO, and 0.71, re- 
spectively. The out-of-plane rigidity is estimate to be 
0.182 eV. These values reveal that the SZOML is much 
softer than graphene. 
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(iii) The SZOML is a wide band gap semiconductor 
with electronic band gap 3.576 eV. This value corre- 
sponds to the first peak in the linear optical spectra 
Ini{ft,ej,j^} — oj. 

We hope these results can be verified by experiments 
in the future. These main results also help us to in- 
fer the possible physical properties of single-walled ZnO 
nanotubes (SWZONTs). Interestingly, there are a lot 
of first-principles investigations on the elastic, electronic, 
and optical properties of SWZONTsii^^i^ although they 
have not been synthesized yet. The geometrical construc- 
tion of SWZONTs from SZOML was discussed in the 
previous workjiSS Here we make some predictions on the 
physical properties of SWZONTs as follows. 

(i) In terms of our experience on the elasticity 
of graphene and carbon nanotubeS )^^'^'^ the in-plane 
Young's modulus, Poisson ratio and out-of-plane bend- 
ing rigidity of SWZONTs should be similar to those of 
SZOML, which is weakly dependent on the chirality of 
the SWZONTs. 

(ii) The piezoelectricity depends on the chirality of 
SWZONTs viewed from the symmetry. Zigzag nanotubes 
possess the largest piezoelectricity while armchair nan- 
otubes have no piezoelectricity. 

(iii) Through zone folding method,-^* we can obtain 
most of phonon dispersion curves of SWZONTs from the 
phonon dispersion relation of SZOML. In particular, the 



sound velocity of twisting mode of SWZONTs is about 
2.801 km/s because this mode corresponds to the in-plane 
tangential acoustic mode of SZOML. 

(iv) Through zone folding method,*^ we can obtain 
the electronic band structure of SWZONTs from the 
band structure of SZOML. Because the SZOML is a 
wide band gap semiconductor, SWZONTs should also 
be wide band gap semiconductors although the band 
gap depends on the their chirality. This point supports 
the result obtained by Elizondo and Mintmire,— but not 
that by Erkog and Kokten.^i Thus electronic property of 
SWZONTs is quite different from that of single- walled 
carbon nanotubes^"* which can be metallic or semicon- 
ducting. 
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